Effect of annealing on the depth profile of hole concentration in (Ga,Mn)As 
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The effect of annealing at 250° C on the carrier depth profile, Mn distribution, electrical conduc- 
tivity, and Curie temperature of (Ga,Mn)As layers with thicknesses > 200 nm, grown by molecular- 
beam epitaxy at low temperatures, is studied by a variety of analytical methods. The vertical 
gradient in hole concentration, revealed by electrochemical capacitance-voltage profiling, is shown 
to play a key role in the understanding of conductivity and magnetization data. The gradient, ba- 
sically already present in as-grown samples, is strongly influenced by post-growth annealing. From 
secondary ion mass spectroscopy it can be concluded that, at least in thick layers, the change in 
carrier depth profile and thus in conductivity is not primarily due to out-diffusion of Mn interstitials 
during annealing. Two alternative possible models are discussed. 
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I. INTRODUCTION 

Semiconductor-based spintronic technology, where 
both the electrical charge and the spin of carriers are 
utilized for signal processing and storage, calls for the 
development of new ferromagnetic materials. The III-V 
dilute magnetic semiconductor (Ga,Mn)As, being com- 
patible with conventional semiconductor technology, is 
considered a potential candidate for spintronics and has 
been intensely studied during the past few years»i*2ii The 
ferromagnetic properties of (Ga,Mn)As, successfully ex- 
plained within the Zener mean-field model^, arise from 
the S = 5/2 spins of Mn atoms incorporated on Ga lattice 
sites. The ferromagnetic Mn-Mn coupling is mediated by 
delocalized or weakly localized holes which are supplied 
by the Mn atoms acting as acceptors on Ga lattice sites. 
The Mn spin system undergoes a ferromagnetic phase 
transition at the Curie temperature Tc which is suggested 
to strongly depend on both, the Mn content and the 
hole concentration. 4 (Ga,Mn)As is grown by molecular- 
beam epitaxy (MBE) at low temperatures (^250°C). In 
most cases, Tq is increased by post-growth annealing at 
temperatures near or even below the growth tempera- 
ture, resulting in values of up to 160 K in (Ga,Mn)As 
single layers, as reported so f ar ^SiL2i2iiSiiiiiS This is 
commonly explained by the removal or rearrangement of 
highly unstable compensating defects such as Mn atoms 
on interstitial lattice sitesiiSaii* 1 ^* 1 ^ 1 ^ The enhancement 
of the ferromagnetism associated with low-temperature 
(LT) annealing is suppressed in the presence of a thin 
GaAs capping layer, indicating that diffusion of defects 
towards the surface plays a crucial role in the anneal- 
ing pro cess It has been revealed by electrochemi- 
cal capacitance- voltage (ECV) profiling and Raman spec- 
troscopy that in general the hole concentration p is not 



constant throughout the (Ga,Mn)As layer but exhibits 
a vertical gradient^ Starting at the GaAs/(Ga,Mn)As 
interface, p monotonously increases and reaches its max- 
imum value near the sample surface. This finding is in 
agreement with the results of spin- wave resonance exper- 
iments which manifest the existence of a gradient in the 
magnetic properties 

For (Ga,Mn)As, it is well known that, within a limited 
temperature range, a slight increase of the growth tem- 
perature results in a pronounced enhancement of both, 
T c and In addition, a significant increase of the sur- 
face temperature, induced by heat radiation from the ef- 
fusion cells, has been detected during LT MBE growth 
of GaAsi 1 ?' 2 ? Therefore, we suppose that besides anneal- 
ing effects during growth a gradual increase of the sur- 
face temperature due do free carrier absorption is likely 
to account for the observed gradient. Even though dis- 
regarded by most other groups so far, the presence of 
the gradient seems to be a general phenomenon, whereas 
its specific profile depends on several parameters such as 
sample thickness and growth conditions. From the p- 
dependence of T c it becomes obvious that a pronounced 
variation in p along the growth direction must have a 
strong impact on the electrical and magnetic properties 
of (GaMn)As. 

In this work, the effect of annealing at 250° C on the 
depth profile of hole concentration in (Ga,Mn)As epilay- 
ers with thicknesses between 0.2 and 1.2 /xm is studied 
by ECV profiling, micro-Raman spectroscopy, conduc- 
tivity measurements, superconducting quantum interfer- 
ence device (SQUID) magnetization measurements, and 
secondary ion mass spectroscopy (SIMS). It is shown that 
the gradient in the carrier density plays a key role in the 
understanding of annealing-induced effects, such as the 
increase in conductivity and Curie temperature. From 
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our SIMS and ECV data we infer that, at least in thick 
layers (~ 1 /an), out-diffusion of Mn interstitials is not 
the dominant mechanism for the annealing-induced en- 
hancement of the hole density, as recently proposed for 
thin layers (< 100 nm)ii2ii£ While most reports in the 
literature are on (Ga,Mn)As epilayers with thicknesses in 
the range of 10-200 nm, where the gradient is probably 
less pronounced, thicker samples, as considered in this 
work, allow to map the depth profiles, e.g., of the hole 
and Mn concentration over larger distances, and there- 
fore to get more detailed information to the understand- 
ing of the defect dynamics. Provided that the fundamen- 
tal diffusion processes during annealing are primarily in- 
dependent of the layer thickness, the findings reported 
in this paper may be helpful in the interpretation of ex- 
perimental data obtained from much thinner samples as 
well. 



II. EXPERIMENTAL DETAILS 

(Ga,Mn)As layers were grown in a PJBER 32 MBE 
machine on In-mounted semi-insulating VGF GaAs(OOl) 
substrates using a conventional Knudsen cell and a hot- 
lip effusion cell to provide the Ga and Mn fluxes, re- 
spectively. A valved arsenic cracker cell was used in 
the non-cracking mode to supply As4 with a maximum 
V/III flux ratio of about 3. First, a GaAs buffer layer 
around 100 nm thick was grown at a temperature of T s 
— 585°C (conventional substrate temperature for GaAs), 
then the growth was interrupted and T s was lowered to 
~250°C. The Mn concentrations in the 0.2-1.2 //m-thick 
(Ga,Mn)As layers under study were determined by flux 
measurements, which have been checked by elastic recoil 
detection measurements (ERD). For details about ERD, 
see Ref . ^ 

ECV analyses were performed using a Bio-Rad PN4200 
profiler. The electrolyte (250 ml aqueous solution of 2.0 
g NaOH + 9.3 g EDTA) is in contact with the semicon- 
ductor forming an electrolyte-semiconductor diode. A 
low-resistance Ohmic contact can be established to the 
(Ga,Mn)As sample without metallization. The poten- 
tial of the sample is measured potentiometrically with 
reference to a saturated calomel electrode. The admit- 
tance Y of the electrolyte-semiconductor contact is de- 
termined by ac measurements under reverse bias at u> c . 
Additionally, the bias voltage V m is wobbled at a consid- 
erably lower modulation frequency to yield the differen- 
tial admittance dY/dV. V m and uj c are selected to ob- 
tain optimum Schottky characteristics of the electrolyte- 
semiconductor contact. The corresponding equivalent 
circuit which is implemented in the analysis software 
consists of the space-charge layer capacitance C and of 
resistances R p and R s connected in parallel and in sc- 
ries, respectively. R s is obtained from the admittance 
data measured at two different frequencies ui c . Using Y, 
dY/dV, and R s , we can determine C and dC/dV yield- 
ing the carrier concentration N(wd) at the edge of the 



depletion region of width Wd according to 
with 

Wd = , (2) 

where e (0.12 nF/m) and e are the permittivity of 
GaAs and the electron charge, respectively. Controlled 
amounts of (Ga,Mn)As are removed in increments of few 
tens of nanometers by passing a dc current I between the 
anodically polarized semiconductor and a carbon counter 
electrode. The removed layer thickness w e is calculated 
from the accumulated transferred charge using Faraday's 
law of electrolysis: 

M { 

where M (144.6), z (6), p (5.36 g/cm 3 ), and F (96490 
Asmol -1 ) denote the molecular weight, effective dissolu- 
tion valence, and density of GaAs, and Faraday's con- 
stant, respectively. Additionally, w e is controlled by me- 
chanical surface tracing. The diode area A of 0.005 to 
0.008 cm 2 is defined by a plastic sealing ring. Normally, A 
is not accurately known at the beginning of an ECV pro- 
filing run, but is routinely measured subsequently. Rel- 
evant recalculation procedures considering the measured 
values of A and the series resistance R s are implemented 
in the original Bio-Rad software of the PN 4200 ECV pro- 
filer. The measured ECV profiles are reproducible within 
an uncertainty in the absolute values of about 15%. More 
information about ECV profiling can be found in Ref.l22l 

Hole concentrations in (Ga,Mn)As can also be esti- 
mated from Raman scattering by coupled plasmon-LO- 
phonon modes ; 23 ' 24 Therefore, micro-Raman measure- 
ments were performed at room temperature (RT) using 
the 514-nm line of an Ar + laser as an excitation source. 
The Raman signals were detected in the backscattering 
configuration z(x, y)z using a DILOR XY 800-mm triple- 
grating spectrometer with a confocal entrance optics and 
a LN 2 -cooled charge-coupled device detector. Further 
experimental details of the micro-Raman measurements 
are given in Ref. |U 

For the electrical measurements, Hall bars with Ti- 
AuPt-Au contacts were prepared on several pieces of the 
cleaved samples without annealing. The contacts were 
checked to be Ohmic with negligibly low resistance. The 
samples were annealed in air using a LINKAM THMS 
600 heating chamber equipped with an electrical feed 
through, which enabled us to perform in situ measure- 
ments of the conductivity. They were mounted on a sil- 
ver block which could be heated electrically or cooled by 
liquid nitrogen over the temperature range from -200 to 
300°C within 2 minutes. 

To measure depth profiles of the Mn fraction, SIMS 
experiments were performed using a commercial Cameca 
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ims4f-E6 spectrometer with Cs + as primary ion beam 
(net impact energy 5.5 keV) at a sputter rate of about 
1 nm/s. In the absence of appropriate calibration stan- 
dards for Mn, the quantitative analysis of the Mn fraction 
refers to the flux measurement of sample B313, which 
gave us a value of 6%. A detailed introduction into the 
SIMS measuring method is given in RefM 

The magnetization measurements were carried out in 
a QUANTUM DESIGN MPMS 5 SQUID magnetometer 
applying an in-plane magnetic field of 5 mT. 



III. RESULTS AND DISCUSSION 

The phenomena discussed in this paper have been qual- 
itatively observed in all (Ga,Mn)As layers grown at V/III 
flux ratios < 3, exhibiting thicknesses between 200 nm 
and 1.2 /im. The specific influence of the V/III flux ra- 
tio on the structural, electric, and magnetic properties of 
(Ga,Mn)As is not subject of this work and will be dis- 
cussed elsewhere. In the following we present experimen- 
tal results obtained from several pieces of a 1.2-^m-thick 
(Ga,Mn)As epilayer with a Mn fraction of 6% (sample 
B313) and of a 240-nm-thick epilayer with a Mn fraction 
of 4.5% (sample B352). These epilayers are representa- 
tive for all other samples investigated so far. 

As mentioned above, ECV profiling reveals the pres- 
ence of a vertical gradient in the hole concentration. Fig- 
urendepicts the ECV profiles of sample B313 before and 
after annealing at 250°C for 30 and 370 min. It is clearly 
seen that the gradient, already present in the as-grown 
sample, is strongly enhanced by post-growth annealing, 
leading to a hole density near the sample surface which is 
almost twice as high as near the GaAs/(Ga,Mn)As inter- 
face. As a consequence, the electric and magnetic prop- 
erties of the (Ga,Mn)As layer are affected in a dramatic 
way, as will be shown below. Whereas the total hole con- 
centration, averaged over the layer thickness, increases 
from (3.8 ± 0.6) x 10 20 cm" 3 in the as-grown sample to 
(4.3 ± 0.6) x 10 20 cm -3 in the 30-min-annealed sample 
and to (4.9 ± 0.7) x 10 20 cm" 3 in the 370-min-annealed 
sample, the local hole density near the sample surface 
almost saturates after annealing for 30 min. The partic- 
ular evolution of the carrier depth profile in the course 
of the annealing seems to corroborate the assumption 
that out-diffusion of compensating defects towards the 
surface accounts for the increase in hole density during 
post-growth annealing>i2iii*i£ 

The ECV profiles of sample B352, recorded before and 
after annealing at 250°C for 30 min, are shown in Fig. El 
The gradient in the as-grown sample is much more pro- 
nounced than in B313 and flattens upon annealing with 
a concomitant increase in the total hole concentration. 
The ECV profiles resemble those of sample B313 in the 
range from the surface down to ~ 240 nm after annealing 
for 30 and 370 min. 

In order to make sure that ECV profiles reflect the 
correct depth distribution of the hole concentration, con- 
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FIG. 1: ECV profiles recorded from sample B313 before and 
after annealing at 250° C for 30 and 370 min. 
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FIG. 2: ECV profiles recorded from sample B352 before and 
after annealing at 250° C for 30 min. Depth-dependent hole 
concentrations determined from Raman spectroscopy are de- 
picted for comparison. 



firmative optical experiments were performed. To this 
end, steplike surface profiles with flat terraces at differ- 
ent depths below the initial surface were prepared on sev- 
eral (Ga,Mn)As layers by wet chemical etching. Micro- 
Raman spectra were taken from each of the individual 
terraces and the carrier densities, obtained from a line- 
shape analysis of the Raman signals, were compared with 
the ECV data. Figure [3 shows, as an example, such a 
comparison for the as-grown sample B352. Within the 
error margin of the two analysis methods, the measured 
depth profiles coincide almost perfectly. 

The Raman spectra corresponding to the data points 
in Fig. [21 are depicted in Fig. [21 The high hole concentra- 
tion in (Ga,Mn) As leads to the formation of a phonon-likc 
coupled mode of the longitudinal optical (LO) phonon 
and the overdamped hole plasmon^& With increasing 
hole concentration, this mode shifts from the frequency 
of the LO phonon to that of the transverse optical (TO) 



4 



1 



g 

;1 



Co 

pi 




J I I I I I I I I I I I LiJ I I I U I I I I I I L 



200 220 240 260 280 300 320 
Raman shift (cm 1 ) 



FIG. 3: Raman spectra recorded from sample B352 for dif- 
ferent etch depths. 



phonon. It is clearly seen in Fig. [31 that the coupled 
mode broadens and shifts to higher frequencies with in- 
creasing etch depth, indicating a decrease in the hole con- 
centration. At an etch depth of 175 nm, the remaining 
layer thickness nearly matches the information depth of 
l/2a ~ 50 nm of the Raman measurement, where a de- 
notes the absorption coefficient at 514 nm wavelength 
of the Ar + laser. Therefore, the narrow Raman line of 
the pure LO-phonon mode in the undoped substrate ap- 
pears at 292 cm -1 . The values for the hole densities were 
obtained from line-shape analyses of the Raman spec- 
tra using a value for the hole mobility of 1.6 cm 2 /Vs. 
The calculated line shapes are drawn as solid lines in 
Fig. 13 Details concerning the calculation of Raman line 
shapes in heavi ly p -doped semiconductors can be found 
in Refs.lHandii 

Under the assumption that the hole mobility does not 
significantly change during the annealing process, an in- 
crease in the hole concentration, as seen from Figs.n an d 
13 should result in an increase in the electrical conductiv- 
ity of the (Ga,Mn)As layer. In fact, such an increase has 
been observed by several authorsi£iL2ii& In the present 
work, the effect of annealing at 250°C on the conductiv- 
ity of sample B313 is shown in Fig.0J 

Monitoring the conductivity by an in situ measurement 
for 400 min reveals that the conductivity monotonously 
rises with a rate that gradually decreases with increasing 
annealing time. After total annealing times of 30, 60, and 
180 min, the annealing process was interrupted for 10 min 
and the temperature was rapidly lowered to 25° C in order 
to probe the RT conductivity (dashed line). During the 
first 30 min, the RT conductivity increases from 91 to 103 
f2 cm . After annealing for 370 min, the RT conduc- 
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FIG. 4: In situ monitoring of the conductivity of sample 
B313 during annealing at 250°C. 



tivity has almost saturated at a value of 116 f2 cm . 
According to these values, the conductivity is enhanced 
by factors of 1.13 and 1.27 upon annealing for 30 and 370 
min, respectively. In contrast to Ref. [6J, no decrease of 
the conductivity for annealing times longer than 2 h is 
observed. Note that the results of our in situ measure- 
ments are in qualitative agreement with those obtained 
from much thinner (Ga,Mn)As epilayers (10-100 nm) at 
lower annealing temperatures (< 200°C)iii& 

We may now compare the annealing-induced increase 
of the conductivity with that in hole concentration. From 
the averaged hole concentrations, derived above from the 
ECV profiles in Fig. ^ we obtain an increase in hole 
density by factors of 1.13 ± 0.3 and 1.29 ± 0.4, which 
are in excellent agreement with the values obtained for 
the conductivity. Thus, the assumption of at most a 
small change in hole mobility during annealing is clearly 
confirmed. From the values of the conductivity and the 
averaged carrier concentrations, effective hole mobilities 
of 1.5 ± 0.3 cm 2 /Vs are deduced for the as-grown as well 
as for the annealed samples. 

Raman spectra, recorded from sample B313 before and 
after annealing at 250° C for 30 and 370 min, are shown 
in Fig. [3] Whereas a strong increase in the hole con- 
centration within the first 30 min of annealing can be 
deduced, the Raman spectrum recorded from the 370- 
min-annealed sample does not significantly differ from 
that of the 30-min-annealed sample. This is in agree- 
ment with the ECV profiles in Fig. ^ which reveal that 
near the sample surface the hole density is only slightly 
enhanced by annealing for more than 30 min. For the 
interpretation of the Raman spectra one should keep in 
mind that the Raman signal, as already mentioned above, 
stems from the near-surface region. The solid lines repre- 
sent model calculations of the Raman line shapes using a 
hole mobility of // = 1.2 cm 2 /Vs and hole concentrations 
of 4.3 x 10 20 cm" 3 and 7.8 x 10 20 cm" 3 for the as-grown 
and the annealed samples, respectively 
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FIG. 5: Raman spectra recorded from sample B313 before 
and after annealing at 250° C. The solid lines are calculated 
line shapes. 

According to the relation^ 

T c oc x x p 1 ^ (4) 

between Curie temperature Tc and hole density p, with 
x denoting the concentration of magnetically active Mn 
ions on Ga sites, an increase in the hole density should 
result in an enhancement of the Curie temperature. In 
fact, such an enhancement is shown by the SQUID mag- 
netization curves in Fig. |HJ recorded from sample B313 
before and after annealing at 250°C for 30 and 370 min. 
Note that the curves are normalized to the values at 5 
K, and thus, cannot be used to obtain information about 
the influence of post-growth annealing on the saturation 
magnetization. We suggest that the extended tails of the 
magnetization curves arise from the vertical gradient of 
the hole density p in the (Ga,Mn)As layer. According 
to Eq. (0J, this gradient results in a depth-dependent 
Curie temperature Tc, and thus the curves in Fig. can 
be viewed as superpositions of individual magnetization 
curves. Then, the values of Tc indicated by arrows have 
to be attributed to the near-surface region, similar to the 
hole densities obtained from the Raman measurements. 
Whereas a Curie temperature Tc of 60 ± 5 K is deduced 
for the as-grown sample, a constant value of 100 ± 5 K 
is obtained for the two annealed samples, yielding an 
increase of Tc by a factor of 1.7 ±0.3. In contrast, an en- 
hancement of Tc by a factor of only 1.2 ± 0.1 would have 
been expected from Eq. and the ECV data measured 
near the sample surface, yielding an increase of the local 
hole density by a factor of 1.7 ±0.5. The discrepancy be- 
tween the two values may be explained by a reduction of 
the number of antiferromagnetically ordered Mn atoms 14 
during annealing. 

The increase in hole density, conductivity, and Curie 
temperature upon post-growth annealing is commonly 
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FIG. 6: Normalized magnetization of sample B313 as a func- 
tion of temperature before and after annealing. 

explained by the removal of compensating defects in 
(Ga,Mn)As. The results of ion channeling experiments 
point to a reduction of Mn interstitials (Mn/), acting as 
compensating double donors* 1 ^ Based on the observation 
of Mn accumulation at the sample surface and on theo- 
retical calculations, this reduction has been traced back 
to an out-diffusion of Mn/ towards the surface followed 
by oxidationpi2*i2i2i In order to verify the latter sugges- 
tion for the thick (Ga,Mn)As layers under study, the Mn 
depth profiles of the as-grown and the annealed samples 
were experimentally determined by SIMS measurements 
recorded from the same sample pieces as used for ECV 
profiling. The Mn profiles are depicted in Figs. [7| and |S| 
for B352 and B313, respectively. Whereas Fig.[7|suggests 
a slight annealing-induced lowering of the Mn fraction 
in the 240-nm-thick sample, no significant difference be- 
tween the three profiles of the 1.2-^m-thick sample can 
be identified in Fig. |H1 Inevitably the question arises if 
an out-diffusion of Mn/ , necessary to completely account 
for the observed ECV profiles, would result in a signifi- 
cant and measureable change in the Mn depth profile at 
all. Therefore, simple considerations are made in the fol- 
lowing to estimate the annealing-induced changes in the 
Mn depth profiles expected for the extreme case that the 
increase in the hole densities in Figs.2] an d0was entirely 
due to the out-diffusion of Mn/ . 

The hole concentration p is given by the density 
[Mnca] of substitutional Mn acceptors (Mnc) minus 
the total density of compensating donors: 

p = [Mn Ga ] -Cx [Mm] -d, (5) 

where [Mn/] denotes the density of Mn/ and d the den- 
sity of all other compensating defects. Note that the 
terms in Eq. (jSJ) are local quantities which may in gen- 
eral vary with depth within the sample. The factor C 
accounts for the fraction of Mn/ acting as donors as well 
as for the corresponding valency (2 for double donors). 
Therefore, the inequality C < 2 holds, where the sign of 
equality applies for the case that all Mn/ act as double 
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out-diffused Mn/ can be caiculated from the conservation 
of the total amount of Mn as follows: 



0.00 0.05 0.10 0.15 0.20 0.25 

Depth z (|j,m) 



FIG. 7: Mn depth profiles of sample B352 before and after 
annealing at 250° C for 30 min measured by SIMS. The dashed 
line represents the calculated Mn profile of the annealed sam- 
ple using Eqs. and JHJ. 
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FIG. 8: Mn depth profiles of sample B313 before and after 
annealing at 250° C for 30 and 370 min measured by SIMS. 
The dashed and dotted lines represent the calculated Mn pro- 
files of the 30-min and 370-min-annealed sample, respectively, 
using Eqs. £J an d ©■ 



donors. The total density [Mn] of Mn atoms is given by 

[Mn] = [Mn Ga ] + [Mm] + [Mn ia ] , (6) 

with [Mriia] denoting the density of electrically inactive 
Mn atoms. Let us now consider the hypothetical case 
that the increase in p, revealed by ECV profiling, is 
solely due to out-diffusion of Mn/ and that inside the 
(Ga,Mn)As layer d, [Mnca], and [Mni a ] remain unaf- 
fected upon annealing. Then, using Eqs. JSJ and ©, the 
local density [Mn] ann of Mn atoms inside the annealed 
sample (z > 0) is given by 



[Mr 



= [Mn] ag - {Pann ~ Pag) / C 



(7) 



where the subscripts ann and ag stand for annealed and 
as-grown, respectively. The increase in the Mn concen- 
tration on the surface (z=0) due to the accumulation of 



([Mn]. 



[Mn] a 



9'z=0 



xS 



dz(p ann -p ag )/C , (8) 



where t denotes the thickness of the as-grown (Ga,Mn)As 
epilayer and S the thickness of the Mn surface layer, 
which is of the order of nm. 

The dashed and dotted lines in Figs. and El represent 
the calculated Mn profiles for the 30-min and 370-min- 
annealed samples, respectively, using Eqs. Q) and © 
with C = 2, and taking into account that a Mn frac- 
tion of 1% corresponds to a concentration of Mn atoms 
of 2.2 x 10 20 cm" 3 . For [Mn] ag , the experimental SIMS 
profiles of the as-grown samples were used, while p ag and 
Pann were taken from Figs.0 an d[21 In the case of sample 
B352 the calculated [Mn] ann curve inside the (Ga,Mn)As 
layer has qualitatively the same form as the measured 
Mn profile, but quantitatively the annealing-induced de- 
crease in Mn fraction is about twice as strong. Note 
that our estimate would yield an even more pronounced 
reduction of the Mn fraction if a value less than 2 had 
been used for the factor C. Thus, it seems that in the 
240-nm-thick sample only a small portion of the enhance- 
ment in hole density is actually due to Mn/ out-diffusion. 
In the case of sample B313 the calculated [Mn] ann curves 
inside the epilayer show a clear deviation from the Mn 
distribution in the as-grown sample, beginning at a depth 
of about 0.5-0.7 iim. and increasing strongly towards the 
surface. This behavior is not even qualitatively reflected 
by the measured Mn profiles. Therefore, we conclude 
that at least in thick (~ 1 /im) samples the observed in- 
crease in hole concentration in the bulk is not primarily 
due to out-diffusion of Mn/. This finding is in agree- 
ment with the low diffusivity of Mn/ at 250°C derived 
in Ref. Hi. Moreover, the calculated narrow peaks at the 
surfaces, arising from the accumulation of out-diffused 
Mn/, are not seen in the measured SIMS profiles. Calcu- 
lation yields {[Mn] ann - [Mn] ag ) z=Q > 6% for B352 and 
([Mn] a nn - [Mn] ag ) z=a > 11% for B313 with S < 10 nm. 
Note however that the SIMS data presented in Fig.[7|and 
[SJmay be too rough to resolve an extremely thin Mn-rich 
oxide layer on the surface, if present. 

Our findings discussed above do not contradict the 
model deduced for much thinner (Ga,Mn)As layers (< 
100 nm) that out-diffusion of Mn/ plays a dominant role. 
In fact, the experimental data available so far suggest 
that, due to the low diffusivity of Mn/ at 250° C, a signif- 
icant out-diffusion of Mn/ occurs in the near-surface re- 
gion, whereas in the bulk the highly unstable Mn/ atoms 
rearrange and form electrically inactive randomly dis- 
tributed precipitates, as proposed in Ref. Within 
the framework of our simple considerations above, this 
means that in the bulk the rise in hole concentration p is 
mainly due to a local increase of [MnjJ and/or [MnGa] 
at the expense of [Mn/] . From the ECV profiles in Fig. 2] 
it follows that this rearrangement of [M n{] does not take 
place homogeneously throughout the whole (Ga,Mn)As 
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layer, leading to an overall upshift of the hole density, 
but exhibits a pronounced dependence on sample depth. 
At present, there is no final explanation for this partic- 
ular evolution of the carrier depth profile on annealing 
which seems to be governed by a diffusion-based mecha- 
nism. Further experiments such as depth-resolved inves- 
tigations on the Mn; distribution and theoretical studies 
have to be performed in the future to clarify this point. 
Yet, we finally close with a brief discussion of two poten- 
tial mechanisms. 

The first mechanism is based on a weak depth- 
dependent diffusion of Mn; towards the surface associ- 
ated possibly with a change of lattice site location. Ini- 
tiated by an out-diffusion of Mn; in the near-surface re- 
gion, Mn; atoms deeper in the bulk successively diffuse 
towards the surface while forming electrically inactive Mn 
clusters and/or MnAs. Due to the low diffusivity of Mn;, 
this process gradually becomes ineffective with increasing 
sample depth. Thus, depending on the annealing time, 
different diffusion-like profiles of Mn; evolve which result 
in the hole density profiles shown in Figs. 2] and [3 This 
explanation also accounts for the fact that, compared to 
the total amount of Mn inside the (Ga,Mn)As layer, only 
a negligibly small fraction of Mn effectively migrates out 
of the bulk. 

The second mechanism is more hypothetical and rests 
upon an assumed out-diffusion of highly mobile, pos- 
sibly compensating, defects other than Mn;, acting at 
least partially as a trigger for the rearrangement of Mn;. 
Arsenic located in interstitial positions (As;) may be a 
potential candidate for such a defect. Similar to LT 
GaAo 28 i 29 i 30 the LT MBE growth of (Ga,Mn)As leads 
to the incorporation of excess As up to 2%, acting at 
least partially as compensating donor defectsAi^ Even 
though controversially discussed, several authors suggest 
that in LT GaAs a considerable fraction of the excess 
As atoms is located in interstitial positions, acting as 
highly mobile defects, while the rest of the excess As is 
believed to be in antisite positions (AsG a )i 31 ' 32 i 33 Asca 
defects are known to remain stable up to 450° C and are 
therefore not expected to play a significant role in the 
physical processes taking place at 250°Ci2i Whereas a 
rearrangement of Mn; could be revealed in (Ga,Mn)As 



by combined channeling Rutherford backscattering and 
by particle-induced x-ray emission experiments 13 , out- 
diffusion of As is difficult to verify. As atoms, which 
diffuse from interstitial sites to the sample surface dur- 
ing annealing, efficiently desorb at the surface and are 
therefore hardly detected by surface sensitive methods. 



IV. SUMMARY 

The depth profile of the hole concentration in thick 
(> 200 nm) MBE-grown (Ga,Mn)As layers, measured by 
ECV profiling, as well as its strong change upon post- 
growth annealing at 250°C has been shown to play a key 
role in the interpretation of conductivity and magnetiza- 
tion data. The annealing-induced increase in the total 
hole concentration, derived from ECV profiling, is in ex- 
cellent quantitative agreement with the change in elec- 
trical conductivity measured in situ during annealing. 
The pronounced enhancement of the hole density near 
the sample surface, confirmed by micro-Raman measure- 
ments, is accompanied by a distinct increase of the Curie 
temperature. The particular evolution of the measured 
ECV profiles under continued annealing suggests that 
diffusion processes play a major role in the post-growth 
annealing of (Ga,Mn)As. From a comparison between 
the ECV profiles and the Mn distributions determined 
by SIMS, it is concluded that, in contrast to thin layers, 
the increase in hole density upon post-growth annealing 
in thick (~ 1 /urn) samples is not primarily due to out- 
diffusion of Mn;. We suppose that a depth-dependent 
rearrangement of Mn; in the bulk, initiated by an out- 
diffusion of Mn; from the near-surface region, accounts 
for the change in the hole-density profile. Alternatively, a 
process based on the out-diffusion of other highly mobile 
defects, such as As;, has been tentatively discussed. 
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